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I .  Introduction 

The  problem  of  determining  lifetimes  for  molecules  in  the  proximity  of 

metal  surfaces  has  been  investigated  for  the  past  twenty  years  due  to  its 

1  2 

fundamental  and  practical  importance.  Both  the  cases  of  a  flat  and  a 
3-5 

structured  surface  have  been  thoroughly  studied,  though  it  seems  to  the 

present  authors  that  for  the  latter  type  of  surfaces,  more  experimental 

* 

studies  are  needed  to  check  the  various  results  obtained  from  theoretical 
2-5 

modeling.  Related  to,  but  in  itself  a  more  interesting  problem,  is  the 
case  where  the  admolecules  are  placed  close  to  a  thin- film  substrate.  Roughly 


speaking,  when  the  thickness  is  thin  enough,  a  thin- film  system  is  no  more 

than  two  surfaces  interferring  with  each  other.  For  the  case  of  a  flat  film, 

2  6 

both  the  problems  of  molecular  fluorescence  lifetimes  ’  and  light-scattering 

7 

by  an  external  source  have  been  studied,  with  the  latter  having  been  extended 

g 

to  optically  nonlinear  films.  A  system  of  perhaps  even  greater  interest  is  a 
corrugated  thin  film,  since  the  two  surface  plasmons  supported  by  the  two 
sides  can  now  couple  with  each  other  leading  to  much  richer  proximity  effects. 
This  cross-coupling  effect  has  been  studied  recently  by  means  of  light¬ 
scattering  experiments  in  which  effects  such  as  the  couplings  into  long-  and 
short-range  surface  plasmons^  ^  and  the  dramatic  enhanced  transparency^ ' 
of  a  metallic  film  have  been  observed.  The  problem  of  molecular  fluorescence 
in  the  vicinity  of  such  a  system  has,  however,  received  less  attention.  Only 
until  very  recently,  being  stimulated  by  a  controversy  among  different  groups 

as  to  whether  cross -coupling  has  really  been  observed  from  the  asymmetry  of 

12 

the  reflectance  dips,  a  group  at  the  University  of  Rochester  has  carried  out 

19-21 

successive  experiments  to  study  the  fluorescence  of  photoresist  molecules 


in  the  vicinity  of  a  Ag  grating  film,  and  it  is  claimed  that  cross-coupling 
has  been  unambiguously  observed  in  their  experiments. 


3 


Motivated  by  the  Rochester  experiments,  we  have  set  forth  here  to 

perform  a  theoretical  investigation  of  the  decay  rates  (inverse  lifetimes)  for 

molecules  near  a  metallic  grating  film  by  extending  our  previous  theory  for  a 

grating  surface.^  We  shall  consider  different  structures  of  the  system, 

including  both  the  asymmetrically  and  symmetrically  (free-standing)  bound 

cases  and  films  with  only  one  surface  roughened,  and  we  shall  make 

9  - 16  19-21 

correlations  with  previous  light-scattering  or  fluorescence  studies 

wherever  possible.  In  the  following,  we  shall  present  in  Sec.  II  a 

phenomonological  theory,  and  numerical  illustrations  of  various  decay-rate 

spectra  will  follow  in  Sec.  III.  In  Sec.  IV  we  shall  discuss  the  possible 

experimental  implications  from  our  calculations. 

Before  we  present  our  formalism,  we  want  to  emphasize  that  whereas  the 

theoretical  treatment  of  the  light- scattering  problem  usually  incorpo*ates  all 

the  dynamics  of  the  problem,  ^  ^  the  case  for  the  calculation  of  the  decay 

rate  problem  in  the  vicinity  of  a  rough  surface  has  been  simplified  to  work  in 

3  4 

the  nonretarded  (static)  limit.  ’  We  have  reevaluated  the  applicability  of 

the  nonretarded  approach  in  terms  of  the  substrate  conductivity  and  the 

molecule-surface  distances,  and  have  emphasized  the  significance  of  applying  a 

dynamic  theory  for  such  decay  rate  calculations.^-  Very  recently,  we  have 

shown  further  that  the  morphology- induced  resonances  in  the  decay  rate 

spectrum  can  only  occur  in  the  dynamic  treatment  but  not  in  the  static 

approach,  due  to  the  neglect  of  radiative  energy  transfer  in  the  static 
22  23 

theory.  ’  Hence,  we  propose  that  to  treat  the  cross-coupling  effects 


(mediated  by  the  grating  period)  in  the  decay-rate  spectrum  for  the 
admolecules,  it  is  necessary  to  impose  a  dynamic  theory  for  the  description, 
which  we  now  proceed  to  develop  in  the  following  section. 
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II.  FORMALISM 

The  geometry  of  our  problem  is  depicted  in  Fig.  1.  where  we  consider  a 
two-level  system  (modeled  by  a  dipole  p)  located  at  a  distance  z  -  d  above  a 
thin  metal  (taken  as  Ag)  film  bounded  by  two  grating  surfaces  located  at  z  -  0 
and  z  -  -t,  respectively.  Thus  the  profile  functions  take  the  simple  forms 


-  r  e 

1  ox 


iQ-jX 


f,  -  e 

2  °2 


iQ2x 


(1) 


respectively,  where  we  shall  assume  small  corrugations  (f  Q.  «  1)  so  that 

5  i 

perturbation  theory  can  be  applied.  The  optical  properties  of  the  three 

media  are  described  by  their  dielectric  functions,  -  £ '  (w)  +  i£”(i*>)  and 

According  to  the  classical  phenomenological  (CP)  approach,  we  must 

calculate  the  reflected  field  (Er)  acting  on  the  emitting  dipole,  which  is 

then  the  cause  of  the  induced  decay,  and  we  then  obtain  the  normalized  decay 

2 

rate  as 


—  -  1  +  — ,  ImG(w) 
70  2k2 


with 
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and  q  being  the  quantum  yield  of  the  emitting  state,  k^  (-  ”)  the  emission 

wave  number  and  the  decay  rate  for  a  free  molecule.  Since  the  details  of 
similar  calculations  for  a  single  grating  surface  have  been  given  in  a 
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previous  article,  here  we  shall  only  sketch  briefly  the  important  steps  in 

obtaining  our  final  results,  where  as  before,  we  shall  consider  the  simple 

case  with  che  dipole  oriented  perpendicular  to  the  film. 

By  generalizing  the  original  integral-equation  theory  for  light- 

scattering  at  a  corrugated  surface,  Mills  and  Maradudin  have  obtained  for 

-*r 

the  corrugated  film  the  /i-th  component  for  the  roughness  contribution  to  E  to 
first-order  perturbation,  given  by^ 


E  (r;u)  - 


16t 


j  d^kj|  e  ^  |  dz’  d^  <kjjW|zz')  l[e^Q(z'+t)  -  e^] 


x  5  (z'  )?1(k|| -kjj0) )  +  [«3-«20(-z'))  5  (z '  +t) 
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v 

« 
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where  ||  denotes  any  vector  on  the  xy-plane,  and  are  Fourier  transforms 

of  Eq.  (1),  0  is  the  Heaviside  step  function,  and  d  is  the  Fourier  transform 

of  the  two-dimensional  "flat  propagator”  obtained  in  Ref.  25.  in  Eq.  (4) 

denotes  the  total  field  for  the  dipole-substrate  problem  in  the  medium  that 

contains  z'  for  the  case  of  a  perfectly  flat  film,  and  expressions  which,  in 

26 

terms  of  the  Green  dyadic,  have  already  been  worked  out  in  the  literature  by 

2  5 

Chance,  Prock  and  Silbey.  As  pointed  out  before,  to  evaluate  the  integral 

of  $(z')  in  Eq.  (4)  where  E^^  may  be  discontinuous  across  the  boundary  of  the 

27 

films,  one  must  adopt  Agarwal's  modifications  and  not  just  take  the  mean 

24 

value  of  the  integrals  at  each  side  of  the  boundary.  It  is  straightforward 

24 

to  show  that  Agarwal’s  modifications  to  the  theory  of  Maradudin  and  Mills 


should  be  applicable  to  multilayered  systems  as  well.  Adopting  further  the 
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and 


7 


4 


x 


wi  ch  k.  -  j  t .  - 

1  1C 


(9) 


Note  that  the  signs  of  and  0 ^  must  be  carefully  assigned  so  as  to  satisfy 

26 

the  boundary  conditions  for  the  fields  at  infinity.  The  functions  and  f£ 
are  defined  as^ 


2Sh. £  0 

f  _  — — LJ - 

1  Sh, £ „  +  h„  £ 


12 


2  1 


(10) 


with 


h^£ -  ih3«2  Can(h2t^ 
h^£2  -  ih2«3  tan(h2t) 
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Hence,  from  Eqs .  (2)  and  (3),  we  finally  obtain  to  first  order  in 
roughness  the  reduced  total  decay  rate  for  the  admolecule ,  given  by  (assuming 


III.  numerical  results 

Here  we  apply  our  general  results  in  Sec.  II  to  several  film  systems  of 
interest,  each  of  which  has  been  studied  previously  in  the  literature  for 
light-scattering  properties. 

A.  Film  with  Only  One  Crating  Surface 

This  system  has  recently  been  studied  by  Friar  and  coworkers^  for  the 

* 

reflectance  from  an  impinging  plane  wave.  To  investigate  our  decay-rate 

problem,  we  simply  set  f  -  0  in  Eq.  (1)  and  are  hence  left  with  only  the 

°1 

term  I0  in  Eq .  (5).  Furthermore,  we  take  -  1.0  and  e  ^  ~  3.6,  and  - 

2  2  8 

£'  +  it"  for  silver  can  be  found  from  various  source  books.  Other 

parameters  are  as  given  in  the  figure  captions . 

Figure  2  shows  the  total  decay  (7,  solid  line)  as  well  as  the  separate 

F  R 

contributions  from  the  flat  (7  ,  dotted  line)  and  rough  (7  ,  dashed  line) 
components  of  the  film.  Here  various  resonance  peaks  can  be  clearly  observed 

corresponding  to  the  single-surface  nonradiative  transfer  (~  3.6  eV) ,  the 

2  22 
flat -film  interference  resonance,  and  the  morphology  induced  resonance  due 

to  the  grating  surface,  respectively. 

Figure  3  shows,  for  the  same  system,  the  total  rate  7  for  different 

thicknesses  of  the  film.  It  is  obvious  that  as  t  -*  «>,  the  single  flat  surface 

result  (dotted  line)  is  recovered,  where  both  the  flat- film  interference  and 

grating- induced  resonance  disappear.  Similar  effects  on  the  light-scattering 

resonances  due  to  the  variation  of  film  thickness  in  this  system  have  been 

observed  in  Ref.  15.  Likewise,  similar  results  are  obtained  if  we  set  f  -  0 

2 

but  i*  0  in  our  expressions  in  Sec.  II,  except  that,  of  course,  in  the 

limit  t  -*  ®,  the  grating  induced  peak  (-  2  eV)  still  shows  up,  in  agreement 

22 

with  our  previous  investigations. 
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B .  Supported  Grating  Film 

F  R 

Figure  4  shows  results  for  7,  7  and  7  ,  where  we  have  assumed  f  —  f 

°1  °2 

-  5  A.  This  asymmetrically -bound  system  -  1,  -  3.6)  is  known  for  its 

"inefficiency"  of  cross-coupling  between  the  two  surface  plasmons  supported  by 

21 

each  grating  surface.  In  fact,  the  resonance  peaks  that  appear  in  Figs.  4 
and  5  re  not  much  different  (except  in  strength)  when  compared  with  the 
results  in  case  A.  This  fact  can  be  easily  understood  in  the  case  of  plane- 
wave  incidence,  where  the  dispersion  relations  for  the  two  surface  plasmons 
are  simpler,  from  which  serious  restrictions  are  found  to  be  imposed  on  the 
coupling  condition  when  the  film  is  bounded  by  two  different  dielectric  media. 

C .  Free-Standing  Grating  Film 

This  is  the  most  efficient  system  for  cross -coupling  studies,  since  \he 

two  surface  plasmons  are  coupled  for  all  energies  in  the  case  of  plane -wave 

21 

incidence.  The  fabrication  of  such  a  film  has  been  described  by  Inagaki  and 
13 

coworkers.  For  simplicity,  we  have  assumed  1  and  f  -  f  -  5k. 

It  is  clearly  seen  from  Fig.  6  that  the  two  surface  plasmons  on  each  grating 
surface  are  coupled  to  the  long-  (u>  -  1.5  eV)  and  short-range  (u>  -  2.0  eV) 
surface  plasmons.  Similar  features  have  been  observed  in  a  light- scattering 


experiment  in  which  the  photoacoustic  technique  has  been  applied  to  monitor 


the  excitation  of  such  cross -coupled  plasmons.  However,  the  splitting  of 

higher-order  peaks  observed  there  does  not  show  up  in  our  calculation  due  to 

the  different  dispersion  relations  which  we  have  here  for  the  coupling  of  the 

dipole  radiations  to  the  film,  instead  of  the  plane-wave  couplings  as  in  the 

light- scattering  case,  where  one  has  multiple  branches  of  dispersion 
13 

relations.  Figure  7  shows  results  for  the  same  system  with  different  film 

thicknesses,  where  again  the  cross-coupling  effect  disappears  for  a  thick 

13 

film,  just  as  in  the  experiments  of  Inagaki  and  coworkers. 
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D.  Remarks  on  Che  Rochester  Experiments 

19-21 

In  three  recent  investigations,  a  group  at  the  University  of 

Rochester  has  performed  beautiful  experiments  in  which  cross-coupling  of  the 
two  surface  plasmons  from  the  sides  of  a  thin  grating  Ag  film  is  unambiguously 
observed,  through  the  scattered  (transmitted)  light  from  the  molecular 
fluorescence  of  a  photoresist  in  the  vicinity  of  the  film.  Both 
asymmetrically*-  and  symmetrically-bound  films  have  been  studied.  We  have 
performed  a  calculation  of  the  decay  rates  [Eqs.  (5)-(19)]  for  one  of  their 
systems  using  their  parameters  for  the  optical  constants  and  the  silver 
grating  of  different  periods,  and  our  results  are  shown  in  Fig.  8.  Although 
we  observe  a  similar  qualitative  feature  to  that  in  Fig.  1  of  Ref.  19,  the 

positions  of  the  peaks  are  quite  different  in  the  two  cases.  We  suspect  that 

v 

the  type  of  experiments  we  are  referring  to  would  be  quite  different '  from  the 

Rochester  experiments,  since  in  the  latter  the  dispersion  relations  are  still 

guided  by  those  derived  from  an  incident  plane  wave,  and  the  fluorescent 

molecules  are  simply  undergoing  Thomson  scattering  by  which  the  incident  light 

from  the  optical  fiber  is  transmitted  through  the  grating  film.  Here,  as 

discussed  above,  we  are  instead  referring  to  the  coupling  between  the  surface 

plasmons  and  the  dipole  fields  of  the  decaying  molecules  in  which  a 

superposition  of  all  the  harmonics  are  involved,  and  hence  we  expect  the 

resonances  to  be  governed  by  different  dispersion  relations.  As  a  matter  of 

21 

fact,  one  of  the  Rochester  experiments  indeed  shows  that  for  the 

s'-  ’metrically-bound  case,  multiple  resonances  due  to  higher-order  dispersions 

do  appear  in  their  transmittance  spectrum,  much  as  in  the  case  of  plane-wave 

13 

scattering  as  studied  by  Inagaki  et  al.  In  fact,  the  peaks  shown  here  for 


the  decay-rate  spectrum  in  Fig.  8  are  not  much  different  from  those  due  to 
single-surface  (thick-film)  couplings  (cf.  Figs.  3  and  5).  Cross -coupling 


12 


into  the  split  modes  of  the  long-  and  short-range  surface  plasmons  can  only  be 
seen  here  if  Che  film  thickness  is  further  decreased,  as  shown  in  Fig.  9. 

We  would  like  to  propose  that  for  the  observation  of  the  peaks 

calculated  in  this  present  work,  fluorescence  experiments  have  to  be  performed 

in  which  the  molecules  are  excited  by  a  short-pulse  laser  and  the  emission 

29 

intensity  is  monitored  by  various  time -resolved  techniques.  In  this 
fashion,  the  decay  rate  values  can  then  be  deduced  by  assuming  appropriate 
forms  for  the  decay  curves  that  fit  the  intensity  data.  In  such  experiments, 
the  pulse  width  of  the  laser  must  be  considerably  shorter  than  the  lifetimes 
(-  I/7)  of  the  admolecules,  and  the  molecule  decaying  from  its  excited  state 
(after  being  "pumped  up"  by  the  laser  pulse)  can  then  be  modeled  as  a 
radiating  dipole.  The  reflected  field  from  the  substrate  film  acting  on  the 
dipole  (which  leads  to  the  induced  decay  of  the  molecule)  can  then  be 
correctly  described  by  the  present  formalism,  and  hence  the  peak  positions 
from  the  present  work  can  be  checked  by  exciting  various  fluorescent  molecules 
of  different  natural  frequencies. 


IV.  CONCLUSION 

We  have  performed  calculations  of  the  decay-rate  spectrum  for 

fluorescent  molecules  near  a  corrugated  metal  film  for  both  supported  and 

unsupported  systems,  and  appropriate  types  of  experiments  have  been  suggested."* 

While  many  qualitative  features  observed  here  are  similar  to  the  case  of 

9  - 16 

light-scattering  analyses  of  such  film  systems,  we  emphasize  that  a 
distinct  feature  of  the  present  problem  is  the  dipole  nature  of  the  radation 
coupled  to  the  surface  plasmons  of  the  film.  Because  of  the  difference 
between  this  coupling  and  the  plane-wave  coupling  to  the  grating  modes,  we 
expect  that  interesting  phenomena  will  arise  if  one  studies  the 
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photoabsorption  of  these  admolecules  in  which  an  external  laser  field  (plane 

wave)  is  introduced  to  drive  the  admolecules  with  both  the  external  and 

molecular  (dipolar)  fields  coupled  to  the  surface  plasmon  modes  at  the  same 

time.  We  have  obtained  some  preliminary  results  for  such  a  problem,  and 

30 

details  will  be  presented  in  a  separate  paper. 
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Figure  Captions 

1.  Configuration  of  the  decay- rate  problem  at  a  corrugated  thin  film. 

2.  Results  for  the  various  reduced  decay  rates  for  a  film  with  one  grating 

surface,  where  we  have  chosen  d  -  200  A,  t  -  100  A,  -  1.0,  <3  -  3.6, 

•  -3  - -1 

f  -  0,  $•  -  50  A  and  Q„  -  10  A  The  solid  line  represents  the 

°1  °2 

total  rate,  and  the  dotted  and  dashed  lines  represent  the  contributions 

m 

from  the  flat  and  rough  components  of  the  film,  respectively. 

3.  Same  as  in  Fig.  2,  but  only  the  total  rates  for  various  film 'thicknesses 
are  shown:  solid  line,  t  -  100  A;  dotted  line  t  -  200  A;  dashed  line,  t  -* 

ao  . 

4.  Results  for  a  supported  ( c ^  —  1.0;  £3  -  3.6)  grating  film,  where  we  have 

-3.-1^ 

taken  d  -  200  A,  t  -  100  A,  f  -  f  -  5  A,  Q,  -  Q.  -  10  A  ,  and  the 

°1  °2  12 
graphical  notations  are  as  given  in  Fig.  2. 

5.  Same  as  in  Fig.  4,  but  only  the  total  rates  for  various  film  thicknesses 

are  shown:  a,  t  -  100  A;  b,  t  -  150  A;  c,  t  -  200  A;  d,  t  -►  ®. 

6.  Results  for  a  free-standing  film  (e^  -  «3  -  1.0),  where  all  other 
parameters  and  notations  are  as  in  Fig.  4. 

7.  Same  as  in  Fig.  6,  but  only  the  total  rates  for  various  film  thicknesses 

are  shown:  solid  line,  t  -  100  A;  dotted  line,  t  -  200  A;  dashed  line, 

t  oo . 

8.  Calculation  of  the  (total)  decay-rate  spectrum  using  the  parameters  given 
in  Ref.  19  with  t  -  500  A.  The  various  graphs  correspond  to  different 
grating  periods  (A  -  2ir/Q) :  a,  A  -  760  run;  b,  A  -  814  nm;  c,  A  -  900  nm; 
d,  A  -  926  nm;  e,  A  -  976  nm.  Note  that  the  logarithmic  scale  is  used  for 
the  reduced  decay  rates . 
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9.  Cross -coupling  into  the  long  and  short-range  surface  plasmons  for  the 

configuration  as  in  Fig.  8  for  thinner  film  thickness  (t  -  150  A).  Only 
the  cases  b  (dashed  line) ,  d  (dotted  line)  and  e  (solid  line)  of  Fig.  8 
are  shown. 
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